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1.1 Introduction 

The last few decades have witnessed a substantial progress in the field of electronic structure of 
materials. Using density functional theory (DFT) |22 25 1 within the local density approxima- 
tion (LDA) or generalized gradient approximation (GGA) [38 1 it is quite a routine to calculate 
the electronic structure of relatively complicated materials containing tens of atoms per unit 
cell. The success of LDA, however, is also accompanied by a number of serious problems. It 
was noticed very early on that when applied to calculate the band structures of s-p semicon- 
ductors and insulators, it was found that the band gaps are systematically underestimated by 
some tens percents. Apart from the too small gaps, the band dispersions are very reasonable. 
This remarkable property of the LDA is still waiting for an explanation since formally there is 
no theoretical justification for identifying the one-particle Kohn-Sham eigenvalues as quasi- 
particle energies observed in photoemission experiments. Applications to alkali metals also 
indicate some problems, albeit less serious. When the band dispersions are compared with 
photoemission data, they are found to be too wide by 10-30 %. Some many -body calcula- 
tions of the electron gas [50 33 1, however, suggest that the band widths are actually widened 
compared with the free-electron values and that the LDA performs better than it is commonly 
believed. If this turned out to be true, photoemission data would presumably need a complete 
revision. In any case, the LDA errors in s-p metals are probably less significant compared to 
the band gap errors in semiconductors and insulators. 

A much more serious problem of the LDA arises when it is applied to calculate the elec- 
tronic structures of so-called "strongly correlated systems". We have to be more precise with 
what we mean by correlations. Even in the electron gas, correlation as conventionally defined 
is rather large. It is as large as exchange that the two almost cancel each other leaving the 
free-electron band essentially unchanged. Thus, it is more appropriate in our case to define 
correlation as anything beyond the LDA rather than anything beyond the Fock exchange since 
the former is usually our starting point in electronic structure calculations of solids. 
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Strongly correlated systems are characterized by partially occupied localized orbitals such 
as found in transition metal oxides or 4f metals. Here the problem is often more of qualita- 
tive rather than quantitative nature. It is often found that the LDA predicts a transition metal 
oxide to be a metal whereas experimentally it is an antiferromagnetic insulator. To cite some 
examples, LaMnC>3, famous for its colossal magnetoresistance, and La2Cu04, a well-known 
parent compound of high-temperature superconductors, are antiferromagnetic insulators but 
predicted to be metals by the LDA |39|. In cases where the LDA does predict the correct 
structure, it is legitimate to ask if the one-particle spectrum is also reproduced correctly. Ac- 
cording to the currently accepted interpretation, transition metal oxides may be classified as 
charge-transfer insulators 1 18 42 1, which are characterized by the presence of occupied and 
unoccupied 3d bands with the oxygen 2p band in between. The gap is then formed by the 
oxygen 2p and unoccupied 3d bands, unlike the gap in LDA, which is formed by the 3d states 
(Mott-Hubbard gap). A more appropriate interpretation is to say that the highest valence state 
is a charge-transfer state: During photoemission a hole is created in the transition metal site 
but due to the strong 3d Coulomb repulsion it is energetically more favourable for the hole 
to hop to the oxygen site despite the cost in energy transfer. A number of experimental data, 
notably 2p core photoemission resonance, suggest that the charge-transfer picture is more ap- 
propriate to describe the electronic structure of transition metal oxides. And of course in the 
case of 4f metals, the LDA, being a one-particle theory, is totally incapable of yielding the 
incoherent part of the spectral function or satellite structures. 

The difficulties encountered by the LDA discussed above have prompted a number of 
attempts at improving the LDA. Notable among these is the GW approximation (GWA), de- 
veloped systematically by Hedin in the early sixties |20|. He showed that the self-energy 
can be formally expanded in powers of the screened interaction W, the lowest term being 
iGW, where G is the Green function. Due to computational difficulties, for a long time the 
applications of the GWA were restricted to the electron gas. With the rapid progress in com- 
puter power, applications to realistic materials eventually became possible about two decades 
ago. Numerous applications to semiconductors and insulators reveal that in most cases the 
GWA 1101 1121 removes a large fraction of the LDA band-gap error. Applications to alkalis 
show band narrowing from the LDA values and account for more than half of the LDA error 
(although controversy about this issue still remains |29|). 

The success of the GWA in sp materials has prompted further applications to more strongly 
correlated systems. For this type of materials the GWA has been found to be less successful. 
For example, GW calculation on nickel 1 8 1 does reproduce the photoemission quasiparticle 
band structure rather well, as compared with the LDA one where the 3d band width is too 
large by about 1 eV, but the too large LDA exchange splitting of 0.6 eV (experimentally 0.3 
eV) remains essentially unchanged. Moreover, the famous 6 eV satellite is not reproduced. 
Application to NiO |9|, a prototype of transition metal oxides, also reveals some shortcom- 
ings. One problem is related to the starting Green's function, usually constructed from the 
LDA Kohn-Sham orbitals and energies. In the LDA the band gap is very small, about 0.2 eV 
compared with the 4 eV experimental band gap. A commonly used procedure of perform- 
ing a one-iteration GW calculation yields about 1 eV gap, much too small. This problem 
is solved by performing a partial self-consistency, where knowledge of the self-energy from 
the previous iteration is used to construct a better starting one-particle Hamiltonian |9|. This 
procedure improves the band gap considerably to a self-consistent value of 5.5 eV and at the 
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same time increases the LDA magnetic moment from 0.9 ^b to about 1.6 /is much closer to 
the experimental value of 1.8 hb ■ However, the GWA maintains the Mott-Hubbard gap, i.e., 
the gap is formed by the 3d states as in the LDA, instead of the charge-transfer gap. In other 
words, the top of the valence band is dominated by the Ni 3d. A more recent calculation using 
a more refined procedure of partial self-consistency has also confirmed these results [ 16 1. The 
problem with the GWA appears to arise from inadequate account of short-range correlations, 
probably not properly treated in the random-phase approximation (RPA). 

Attempts at improving the LDA to treat strongly correlated systems were initiated by the 
LDA+U method 0B]|3O]|3, which introduces on top of the LDA Hamiltonian a Hubbard U 
term and a double-counting correction term, usually applied to partially filled 3d or 4f shells. 
The LDA+U method is essentially a Hartree-Fock approximation to the LDA+U Hamitonian. 
In the LDA, the Kohn-Sham potential does not explicitly distinguish between occupied and 
unoccupied orbitals so that they experience the same potential. In, for example, transition 
metal oxides, where the 3d orbitals are partially occupied, this leads to metallicity or to under- 
estimation of the band gap. The LDA+U cures this problem by approximately pushing down 
the occupied orbitals by U/2 and pushing up the unoccupied orbitals by U/2, creating a lower 
and upper Hubbard band, thus opening up a gap of the order of the Hubbard U. The LDA+U 
method has been successfully applied to late transition metal oxides, rare earth compounds 
such as CeSb, as well as to problems involving metal-insulator transition and charge-orbital 
ordering. 

More recently, the idea of the LDA+U was extended further by treating the Hubbard U 
term in a more sophisticated fashion utilizing the dynamical mean-field theory (DMFT) 1 19 1. 
The DMFT is remarkably well suited for treating systems with strong on-site correlations be- 
cause the on-site electronic Coulomb interactions are summed to all orders. This is achieved 
by using a mapping onto a self-consistent quantum impurity problem, thereby including the 
effects of the surrounding in a mean-field approximation. The strength of the DMFT is its 
ability to properly describe Mott phenomenon or the formation of local moments, which is 
the key to understanding many physical properties in strongly correlated materials. The com- 
bination of LDA and DMFT takes advantage of the first-principles nature of LDA while at 
the same time incorporates local correlation effects not properly treated within the LDA. The 
LDA+DMFT method [7 6 31 1 has been successfully applied to a number of systems by now. 

In both the LDA+U and LDA+DMFT methods, two fundamental problems remain un- 
addressed. First, the Hubbard U is usually treated as a parameter, and second, the Hubbard 
U term contains interaction already included in the LDA but it is not clear how to take into 
account this double-counting term in a precise way. Thus, a truly first-principles theory for 
treating strongly correlated systems is still lacking. In this article, we describe a dynamical 
mean-field approach for calculating the electronic structure of strongly correlated materials 
from first-principles 1 1 31 1461 . The DMFT is combined with the GW method, which enables 
one to treat strong interaction effects |28|. One of the main features of the new scheme is 
that the Hubbard U is calculated from first principles through a self-consistency requirement 
on the on-site screened Coulomb interaction, analogous to the self-consistency in the local 
Green's function in the DMFT. Since the GWA has an explicit diagrammatic representation, 
the on-site contribution of the GW self-energy can be readily identified and the scheme then 
allows for a precise double-counting correction. 

In the next two sections, we will give a summary of the GWA and DMFT, describing 
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their main features. In the fourth section we lay out the GW+DMFT scheme, followed by a 
simplified application of the scheme to the excitation spectrum of nickel. Finally we discuss 
some future challenges and directions. 

1.2 The GW Approximation 
1.2.1 Theory 

It can be shown that the self-energy may be expressed as l20l 

£(1,2) = -* J d3dAv(l,4)G(l,3) 6G ~^ 2) (1.1) 

where v is the bare Coulomb interaction, G is the Green function and <f> is an external time- 
dependent probing field. We have used the short-hand notation 1 = (xiti). From the equation 
of motion of the Green function 

G- 1 = & - H Q - £ (1.2) 
at 



H Q = h Q + <j> + V H (1.3) 
ho is the kinetic energy and Vjj is the Hartree potential. We then obtain 
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where e 1 is the inverse dielectric matrix. The GWA is obtained by neglecting the vertex 
correction <5£/(50, which is the last term in HAl . This is just the random-phase approximation 
(RPA) for e _1 . This leads to 



£(1, 2) = iG(l, 2)^(1,2) (1.5) 
where we have defined the screened Coulomb interaction W by 

W(l,2) = J d3v(l,3)e- 1 (3, 2) (1.6) 
The RPA dielectric function is given by 



e=l-vP 



(1.7) 
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with the Green function constructed from a one-particle band structure {i/ji, £j}. The factor 
of 2 arises from the sum over spin variables. In frequency space, the self-energy in the GWA 
takes the form 



We have so far described the zero temperature formalism. For finite temperature we have 



In the Green function language, the Fock exchange operator in the Hartree-Fock approxi- 
mation (HFA) can be written as iGv. We may therefore regard the GWA as a generalization 
of the HFA, where the bare Coulomb interaction v is replaced by a screened interaction W. 
We may also think of the GWA as a mapping to a polaron problem where the electrons are 
coupled to some bosonic excitations (e.g., plasmons) and the parameters in this model are 
obtained from first-principles calculations. 

The replacement of v by W is an important step in solids where screening effects are 
generally rather large relative to exchange, especially in metals. For example, in the electron 
gas, within the GWA exchange and correlation are approximately equal in magnitude, to a 
large extent cancelling each other, modifying the free-electron dispersion slightly. But also 
in molecules, accurate calculations of the excitation spectrum cannot neglect the effects of 
correlations or screening. The GWA is physically sound because it is qualitatively correct in 
some limiting cases 1211 . 

1.2.2 The GW approximation in practice 

The quality of the GWA may be seen in Figure dl.lt . where a plot of band gaps of a number of 
well known semiconductors and insulators is displayed. It is clear from the plot that the LDA 
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Figure 1.1: Band gaps of some selected semiconductors and insulators calculated within the GWA 
compared with the LDA and experimental values. The GW data are taken from 1 26 1. 



systematically underestimates the band gaps and that the GWA substantially improves the 
LDA band gaps. It has been found that for some materials, like MgO and InN, significant error 
still remains within the GWA. The reason for the discrepancy has not been understood well. 
One possible explanation is that the result of the one-iteration GW calculation may depend 
on the starting one-particle band structure. For example, in the case of InN, the starting LDA 
band structure has no gap. This may produce a metal-like (over)screened interaction W which 
fails to open up a gap or yields too small gap in the GW calculation. A similar behaviour is 
also found in the more extreme case of NiO, where a one-iteration GW calculation only yields 
a gap of about 1 eV starting from an LDA gap of 0.2 eV (the experimental gap is 4 eV) I9l ll0l . 

The problems with the GWA arise when it is applied to strongly correlated systems. Ap- 
plication to ferromagnetic nickel 1 8 1 illustrates some of the difficulties with the GWA. Starting 
from the LDA band structure, a one-iteration GW calculation does improve significantly the 
LDA band structure. In particular it reduces the too large 3d band width bringing it into much 
better agreement with photoemission data. However, the too large LDA exchange splitting 
(0.6 eV compared with experimental value of 0.3 eV) remains essentially unchanged. More- 
over, the famous 6 eV satellite, which is of course missing in the LDA, is not reproduced. 
These problems point to deficiencies in the GWA in describing short-range correlations since 
we expect that both exchange splitting and satellite structure are influenced by on-site inter- 
actions. In the case of exchange splitting, long-range screening also plays a role in reducing 
the HF value and the problem with the exchange splitting indicates a lack of spin-dependent 
interaction in the GWA: In the GWA the spin dependence only enters in G but not in W . 

Application to NiO, the prototype of Mott-Hubbard transition metal oxides, reveals an- 
other difficulty with the one-iteration GWA. As already mentioned previously, when the start- 
ing band structure is far from the experimental quasiparticle band structure, a one-iteration 
GW calculation may not be sufficient. This problem may be circumvented by performing 
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a partial self-consistent calculation in which the self-energy from the previous iteration at a 
given energy, such as the Fermi energy of the centre of the band of interest, is used to con- 
struct a new set of one-particle orbitals. This procedure is continued to self-consistency such 
that the starting one -particle band structure gives zero self-energy correction l9l llQ|[T6l . A 
more serious problem, however, is describing the charge-transfer character of the top of the 
valence band. The GWA essentially still maintains the Mott-Hubbard band gap as in the LDA, 
i.e., the top of the valence band is mainly of 3d character rather than the charge-transfer char- 
acter dominated by the 2p oxygen hole. As in nickel, the problem with the satellite arises 
again. Depending on the starting band structure, a satellite may be reproduced albeit at a too 
high energy. Thus there is a strong need for improving the short-range correlations in the 
GWA which may be achieved by using a suitable approach based on the dynamical mean-field 
theory described in the next section. 



Dynamical mean field theory (DMFT) |19| has originally been developed within the con- 
text of models for correlated fermions on a lattice where it has proven very successful for 
determining the phase diagrams or for calculations of excited states properties. It is a non- 
perturbative method and as such appropriate for systems with any strength of the interaction. 
In recent years, combinations of DMFT with band structure theory, in particular Density func- 
tional theory with the local density approximation (LDA) have emerged. The idea is to correct 
for shortcomings of DFT-LDA due to strong Coulomb interactions and localization (or par- 
tial localization) phenomena that cause effects very different from a homogeneous itinerant 
behaviour. Such signatures of correlations are well-known in transition metal oxides or f- 
electron systems but are also present in several elemental transition metals. 

The application of DMFT to real solids relies on a representability conjecture assuming 
that local quantities such as for example the local Green's function or self-energy of a solid 
can be calculated from a local impurity model, that is one can find a dynamical mean field 
Qo and a Hubbard parameter U, such that the Green's function calculated from the effective 
action 



coincides with the local Green's function of the solid. For a model of correlated fermions on 
a lattice with infinite coordination number this conjecture can rigorously be proven: it is a 
consequence of the absence of any non-local contributions to the self-energy of the system. 
For a real solid the situation is somewhat more complicated, not only due to the finite co- 
ordination number but also to the difficulty to define the notion of locality. This notion is 
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crucial at both stages, for the construction of the impurity model, where long-range Coulomb 
interactions are mimicked by local Hubbard parameters 2 , and for the resolution of the model 
within DMFT, which approximates the full self-energy of the model by a local quantity. Ap- 
plications of DMFT to electronic structure calculations (e.g. the LDA+DMFT method) are 
therefore always defined within a specific basis set using localized basis functions. Within an 
LMTO implementation for example locality can naturally be defined as referring to the same 
muffin tin sphere. This amounts to defining matrix elements Glr.l'R' (iw) of the full Green's 
function 

G(r,r',iw) = xl R (r)G L R,L'R'(HxL'R'( r ') 

LL'RR' 

and assuming that its local, that is "on-sphere" part equals the Green's function of the local 
impurity model (11.1 3i . Here R, R' denote the coordinates of the centres of the muffin tin 
spheres, while r, r' can take any values. The index L = (n,l,m) regroups all radial and 
angular quantum numbers. The dynamical mean field Qo in Jl . 1 3i has to be determined in such 
a way that the Green's function Gi mpur it y L,L',a of the impurity model Ea. dl. 131 coincides 
with Glr.l'R' (ioj) if the impurity model self-energy is used as an estimate for the true self- 
energy of the solid. This self-consistency condition reads 

GimpurityiiUn) = 2_j + A* _ H o(k) - S(iw„)) _1 

k 

where E, Ho and G are matrices in orbital and spin space, and iu> + /i is a matrix proportional 
to the unit matrix in that space. 

Together with Jl . 1 3i this defines the DMFT equations that have to be solved self-consistently. 
Note that the main approximation of DMFT is hidden in the self-consistency condition where 
the local self-energy has been promoted to the full lattice self-energy. 

The representability assumption can actually be extended to other quantities of a solid than 
its local Green's function and self-energy. In "extended DMFT" 1441 1271 l24l 1431 e.g. a two 
particle correlation function is calculated and can then be used in order to represent the local 
screened Coulomb interaction W of the solid. This is the starting point of the "GW+DMFT" 
scheme described in section 6. 



1.3.1 DMFT in practice 

Combinations of DFT-LDA and DMFT, so-called "LDA+DMFT" techniques Q, have so 
far been applied to transition metals (Fe, Ni, Mn) and their oxides (e.g. La/YTi03, V2O3, 
Sr/CaVOs) as well as elemental f-electron materials (Pu, Ce) and their compounds. In the 



2 For a discussion of the appropriateness of local Hubbard parameters see II II . 
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most general formulation, one starts from a many-body Hamiltonian of the form 

H= ^ {H-im,Vm> - H d °)<4mv a i'm'cr (1-14) 
{wicr} 

imm'cr 



2 . 



where H is the effective Kohn-Sham-Hamiltonian derived from a self-consistent DFT- 
LDA calculation. This one-particle Hamiltonian is then corrected by Hubbard terms for direct 
and exchange interactions for the "correlated" orbitals, e.g. d or / orbitals. In order to avoid 
double counting of the Coulomb interactions for these orbitals, a correction term H dc is sub- 
tracted from the original LDA Hamiltonian. The resulting Hamiltonian d 1 . 1 4i is then treated 
within dynamical mean field theory by assuming that the many-body self-energy associated 
with the Hubbard interaction terms can be calculated from a multi-band impurity model. 

This general scheme can be simplified in specific cases, e.g. in systems with a separation of 
the correlated bands from the "uncorrelated" ones, an effective model of the correlated bands 
can be constructed; symmetries of the crystal structure can be used to reduce the number of 
components of the self-energy etc. 

Despite the huge progress made in the understanding of the electronic structure of cor- 
related materials thanks to such LDA+DMFT schemes, certain conceptual problems remain 
open: These are related to the choice of the Hubbard interaction parameters and to the dou- 
ble counting corrections. An a priori choice of which orbitals are treated as correlated and 
which orbitals are left uncorrelated has to be made, and the values of U and J have to be fixed. 
Attempts of calculating these parameters from constrained LDA techniques are appealing in 
the sense that one can avoid introducing external parameters to the theory, but suffer from 
the conceptual drawback in that screening is taken into account in a static manner only II II . 
Finally, the double counting terms are necessarily ill defined due to the impossibility to single 
out in the LDA treatment contributions to the interactions stemming from specific orbitals. 
These drawbacks of LDA+DMFT provide a strong motivation to attempt the construction 
of an electronic structure method for correlated materials beyond combinations of LDA and 
DMFT 



1.4 GW+DMFT 

The idea behind combining the GWA and the DMFT is to take advantage of the strong features 
of the two theories. The GWA, being based on RPA screening, is capable of taking into 
account long-range correlations but does not describe properly short-range correlations. On 
the other hand, the strength of the DMFT is its ability to describe on-site correlations by 
means of a mapping of the many-body problem onto an impurity problem where the on- 
site interactions are summed to all orders. While the GWA is a fully first-principles theory, 
the DMFT has been traditionally used in conjunction with a Hubbard model. Clearly, it is 
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desirable to combine the two theories into a consistent theory where the parameters in the 
Hubbard model are determined self-consistently from first-principles via the GWA. 

Recently, first steps have been undertaken towards a combination of the GWA and DMFT, 
both in a model context |46| and within the framework of realistic electronic structure calcu- 
lations 1131 . The basic physical idea of GW+DMFT is to separate the lattice into an on-site 
part and the rest. The on-site self-energy is taken to be the impurity self-energy calculated 
by the DMFT and the off-site self-energy is calculated by the GWA. Viewed from the GWA, 
we replace the on-site GW self-energy by that of DMFT, correcting the GW treatment of on- 
site correlations. Viewed from the DMFT, we add off-site contributions to the self-energy 
approximated within the GWA, giving a momentum dependent self-energy. 

The impurity problem contains a Hubbard interaction U that is usually treated as a pa- 
rameter. In order to calculate the U, we introduce a self-consistency condition that the U 
screened by the effective bath in the impurity model be equal to the local projection of the 
global screened interaction W. This condition is complementary to the condition imposed in 
the DMFT that the impurity Green's function be equal to the local Green's function. 

The above physical ideas can be formulated in a precise way using the free-energy func- 
tional of Luttinger and Ward (LW). A generalization of the original LW functional takes the 
form 1 1 1 (see also 1 15 1) 



r(G, W) ^Trh^G-Tr[(G H 1 -G- 1 )G]-^TrlnW+^Tr[(v- 1 -W- 1 )W]+'i'[G,W} 

(1.15) 

Gft 1 = iLo n + /i + V 2 /2 — Vh corresponds to the Hartree Green's function with Vh being 
the Hartree potential. The functional ^[G, W] is a generalization of the original LW $[G] 
functional, whose derivative with respect to G gives the self-energy. A more general derivation 
of J 1 . 1 51 using a Hubbard-Stratonovich transformation and a Legendre transformation with 
respect to both G and W may be found in a later work 1151 . It is straightforward to verify that 
at equilibrium the stationarity of T yields 

£-0^-.--P,P~2» 
The functional is divided into on-site and off-site components: 

* = ^>°Jl' sl te [G RR \ W RR '} + *Zp SUe [ GRR , wRR \ (1-18) 
where R denotes a lattice site. In the GWA the functional is given by 



$gw[G, W] = -GWG 



(1.19) 
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The impurity part \&, mp is generated from a local quantum impurity problem defined on a 
single atomic site with an effective action 

S= J drdr' [-J2ci(r)g^ L ,(r-r')c L ,(r') 

+ \Y.'- c U t ) c ^) ■U LiL2LsU {t-t') :c+ (t') C l 4 (t) : (1.20) 

The double dots denote normal ordering and L refers to an orbital of angular momentum L 
on a given sphere where the impurity problem is defined. These orbitals are usually partially 
filled localized 3d or 4f orbitals. 

The GW+DMFT set of equations can now be readily derived from dl.l5> . dl.l8> . and 
J 1 . 1 91 by taking functional derivatives of \& with respect to G and W as in dl.l6t and dl . 17l >: 

E = Eg£(l - Sjw) + Z™5 RW (1.21) 



P = Pgw (1 - 8rr>) + P£? p 6rr, (1.22) 

In practice the self-energy is expanded in some basis set {4>l} localized in a site. The polar- 
ization function on the other hand is expanded in a set of two-particle basis functions {<fiL<pL'} 
(product basis) since the polarization corresponds to a two-particle propagator. For example, 
when using the linear muffin-tin orbital (LMTO) band-structure method, the product basis 
consists of products of LMTO's. These product functions are generally linearly dependent 
and a new set of optimized product basis (OPB) [ 10 1 is constructed by forming linear combi- 
nations of product functions, eliminating the linear dependencies. We denote the OPB set by 
B a = 4>l<I>l'c'll'- T° summarize, one-particle quantities like G and E are expanded in 

{4>l} whereas two-particle quantities such as P and W are expanded in the OPB set {B a }. It 
is important to note that the number of {B a } is generally smaller than the number of {4>l 4>l' } 
so that quantities expressed in {B a } can be expressed in {4>l4>l' }, but not vice versa. In mo- 
mentum space, equations I1.2H and dl.22t read 



E Li (k, icu n ) = E^V (k, *w« ) " 22 ^gw (k, icon ) + (tw n ) (1.23) 



P a "(k,lU; n ) = J^(k,iw„) ~J2 P Gw(^^n)+P^ p (^n) (1.24) 

k 

The second terms in the above two equations remove the on-site contributions of the GW 
self-energy and polarization, which are already included in Ei mp and Pi mp . 

We are now in a position to outline the self-consistency loop which determines Q and U 
as well as the full G and W self-consistently. 
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• The impurity problem dl.20t is solved, for a given choice of Weiss field Qll' and Hub- 
bard interaction U a p: the "impurity" Green's function 

G^ p = -{T T c L {r)c+{r')) s (1.25) 
is calculated, together with the impurity self-energy 

Sj mp = Styi mp /SGi m p = Q^ 1 - G^ mp . (1-26) 

The two-particle correlation function 

Xl ± l 2 l 3 l 4 = (■ c ] Li {t)c L2 {t) :: c[ 3 (t')c L4 (V ) :} s (1.27) 
must also be evaluated. 

• The impurity effective interaction is constructed as follows: 

W£? p = Z4/3- £ Y. U ^°^L 2 XL 1 L 2 L 3 L i [Ol iLi rU S p (1.28) 
L\---L& 75 

Here all quantities are evaluated at the same frequency 3 and Oj^ L is the overlap matrix 
(-B 7 1 0l The polarization operator of the impurity problem is then obtained as: 

P imp = -2S% mp /SW imp = U- 1 - Wrf p , (1.29) 

where the matrix inversions are performed in the OPB {B a }. 

• From Eqs. dl.23> and l l 1.241 the full k-dependent Green's function G(k, iuo n ) and ef- 
fective interaction W(q, iv n ) can be constructed. The self-consistency condition is ob- 
tained, as in the usual DMFT context, by requiring that the on-site components of these 
quantities coincide with Gi mp and Wi mp . In practice, this is done by computing the 
on-site quantities 

G loc (iu n ) = 5^[GH _1 (k,»Wn) - ^(K^n)]- 1 (1.30) 

k 

Wi oc (iu n ) = ^[Vq- 1 - P(q, iVnT 1 (1.31) 
q 

and using them to update the Weiss dynamical mean field Q and the impurity model 
interaction U according to: 

Q~ = Gj oc + S imp (1-32) 
U ~ X = W^l + P imp (1.33) 



3 Note that \L i... L 4 does not denote the matrix element < L1L2 >> but is rather defined by x( r > r ') 

T. Ll ..L 4 <t>l 1 l r )<t>t 2 (r)XLi...Li4>L 3 {r')<t>L 4 (r'). 
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This cycle is iterated until self-consistency for Q and U is obtained (as well as on G, W, 
Yj xc and P). When self-consistency is reached, Gi mp — Gi oc and Wi mp = Wi oc . This implies 
that at self-consistency, the second term in dl.231 can be rewritten as (in imaginary-time) 

££^(k,r) = - W^ L ' L hr)G^{r) (1.34) 

k L 1 L' 1 

This shows that the local or on-site contribution of the GW self-energy is precisely sub- 
tracted out, thus avoiding double counting. Eventually, self-consistency over the local elec- 
tronic density can also be implemented, (in a similar way as in LDA+DMFT I40l H-Il t by 
recalculating p(r) from the Green's function at the end of the convergence cycle above, and 
constructing an updated Hartree potential. This new density is used as an input of a new GW 
calculation, and convergence over this external loop must be reached. While implementing 
self-consistency within the GWA is known to yield unsatisfactory spectra 1231 . we expect a 
more favourable situation in the proposed GW+DMFT scheme since part of the interaction 
effects are treated to all orders. 



1.4.1 Simplified implementation of GW+DMFT and application to 
ferromagnetic nickel 

The full implementation of the proposed approach in a fully dynamical and self-consistent 
manner is at the present stage computationally very demanding and we regard it as a major 
challenge for future research. Here, we apply a simplified scheme of the approach 1 13 1 to the 
electronic structure of nickel in order to demonstrate its feasibility and potential. The main 
simplifications made are: 

1 . The DMFT local treatment is applied only to the d-orbitals, and we replace the dynamical 
impurity problem by its static limit, solving the impurity model dl.20> for a frequency- 
independent U = U(uj = 0). 

2. We perform a one-iteration GW calculation in the form 1 10 1 : Sgw = Glda'W[Glda], 
from which the off-site part of the self-energy is obtained. 

The local Green's function is taken to be 

Gf oc (icOn) = Y^iGH-'ik,^) - Z°J£- site (1.35) 

k 

~ {^imp,<T ~ ~^Tr a Y]i m p ta (0) + Vxc site )} 

Thus, the off-site part is obtained from the GW self-energy whereas the on-site part is derived 
from the impurity self-energy with a double-counting correction of the form proposed in 1 32 1. 

We have performed finite temperature GW and LDA+DMFT calculations (within the 
LMTO-ASA|2| with 29 irreducible k-points) for ferromagnetic nickel (lattice constant 6.654 
a.u.), using 4s4p3d4f states, at the Matsubara frequencies ioj n corresponding to T = 630A", 
just below the Curie temperature. The GW self-energy is calculated from a paramagnetic 
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Figure 1.2: Partial density of states pf 3d orbitals of nickel (solid/dashed lines give the majority/minority 
spin contribution) as obtained from the combination of GW and DMFT 1131 . For comparison with LDA 
and LDA+DMFT results, see |32|; for experimental spectra, see 1351 . 



Green's function, leaving the spin-dependence to the impurity self-energy. The resulting self- 
energies are inserted into Eq. (11.351 . which is then used to calculate a new Weiss field ac- 
cording to dl .32i . The Green's function G^ c (r) is recalculated from the impurity effective 
action by QMC and analytically continued using the Maximum Entropy algorithm. The re- 
sulting spectral function is plotted in Fig.l ll.2l l. Comparison with the LDA+DMFT results in 
1 32 1 shows that the good description of the satellite structure, exchange splitting and band 
narrowing is indeed retained within the (simplified) GW+DMFT scheme. We have also calcu- 
lated the quasiparticle band structure, from the poles of dl .351 . after linearization of S(k, iu> n ) 
around the Fermi level 4 . Fig. dl .31 shows a comparison of GW+DMFT with the LDA and 
experimental band structure. It is seen that GW+DMFT correctly yields the bandwidth reduc- 
tion compared to the (too large) LDA value and renormalizes the bands in a (k-dependent) 
manner. 

Because of the static approximation 3), we could not implement self-consistency on Wi oc 
(Eq. I1.3H ). We chose the value of U{u> = 0) (~ 3.2eV) by calculating the correlation 
function x and ensuring that Eq. dl .28b is fulfilled at u> = 0, given the GW value for Wi oc (oj = 
0) (~ 2.2eV for Nickel l45l l. 

1.5 Conclusions 

The proposed GW+DMFT scheme avoids the conceptual problems inherent to "LDA+DMFT" 
methods, such as double counting corrections and the use of Hubbard parameters assigned to 
correlated orbitals. The notion of self-consistency condition on the on-site Green's function in 
the DMFT is extended to the screened interaction. Analogous to the usual condition that the 



4 Note however that this linearization is no longer meaningful at energies far away from the Fermi level. We 
therefore use the unrenormalized value for the quasi-particle residue for the s-band (Z s = 1). 
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Figure 1.3: Band structure of nickel (majority and minority spins) from GW+DMFT scheme 1 13 1 (cir- 
cles) in comparison to the LDA band structure (dashed lines) and experiments 1 14 1 (triangles down) and 
1351 (triangles up). 



impurity Green's function be equal to the on-site Green's function, we demand that the on-site 
screened Hubbard U of the impurity be equal to the on-site projection of the global screened 
interaction W. In this fashion, the Hubbard U is determined from first-principles. Since the 
GWA has a well-defined diagrammatic interpretation, it is also possible to precisely take into 
account the double-counting correction. 

A number of issues are of immediate importance. Solving impurity models with frequency- 
dependent interaction parameters 1 34 17 46 1 is one of the most urgent tasks as well as studies 
of various possible self-consistency schemes. To study these aspects, we are now carrying out 
GW+DMFT calculations on the electron gas. Applications to real materials are both theoreti- 
cally and computationally very challenging. Extension to the two-particle Green's function is 
another field for future research [371. 
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